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Since the first quantum ghost imaging (QGI) experiment in 1995, many QGI schemes have been 
put forward. However, the position-position or momentum-momentum correlation required in these 
QGI schemes cannot be distributed over optical fibers, which limits their large geographical applica¬ 
tions. In this paper, we propose and demonstrate a scheme for long distance QGI utilizing frequency 
correlated photon pairs. In this scheme, the frequency correlation is transformed to the correlation 
between the illuminating position of one photon and the arrival time of the other photon, by which 
QGI can be realized in the time domain. Since frequency correlation can be preserved when the 
photon pairs are distributed over optical fibers, this scheme provides a way to realize long-distance 
QGI over large geographical scale. In the experiment, long distance QGI over 50 km optical fibers 
has been demonstrated. 


Quantum ghost imaging (QGI) [1] has attracted much 
attention in last two decades due to its abundance in 
physics and potential for quantum communication and 
quantum sensing. Originally, QGI was realized by photon 
pairs generated by spontaneous parametric down con¬ 
version (SPDC) in nonlinear crystals. The momentum- 
momentum correlation in the photon pairs was utilized in 
QGI to realize imaging in a nonlocal manner. After that, 
many ghost imaging schemes have been proposed [2-10]. 
The concept of the ghost imaging has been deeply devel¬ 
oped “from quantum to classical to computational” [11]. 
While, the motivation of researches on ghost imaging 
also developed from the interest on fundamental physics 
[12,13] to applications such as information transmission 
and sensing [14]. It has been reported that ghost imaging 
can be used in optical information encryption and trans¬ 
mission [14]. However, it is still a problem how to realize 
QGI over large geographical scale by techniques compat¬ 
ible with optical communication and networks, since sin¬ 
gle mode fiber (SMF) used in optical communication does 
not support the distribution of the position-position or 
momentum-momentum correlation required in previous 
QGI schemes. 

In this paper, we propose and demonstrate a scheme 
realizing long-distance QGI over optical fibers, based on 
frequency correlated photon pairs. It’s well known that 
the frequency is a stable degree of freedom (DOF) of pho¬ 
tons traveling in optical fibers. The frequency correlation 
can be preserved when the photon pairs are distributed 
over optical fibers. After long distance distribution over 
optical fibers, the two photons in the pair are sent to Al¬ 
ice and Bob sides, respectively. At Alice side, the photons 
are spatially dispersed to different directions according to 
their frequencies by a spatial dispersion component, such 
as a grating, then illuminate the object on different posi¬ 
tions along a line. At Bob side, a component with large 
temporal dispersion is used to change the arrival time 
of idler photons with different frequencies when detected 
by a single photon detector (SPD). Hence, the frequency 
correlation in the initial photon pairs is transformed to 
the correlation of the illuminating positions of the pho¬ 


tons at Alice side and the arrival time of the photons at 
Bob Side. Based on this correlation, the image along the 
illuminating line on the object can be obtained in time 
domain by the coincidence measurement, realizing one¬ 
dimensional QGI. Two-dimensional imaging can also be 
realized by step-moving the object, realizing a function of 
long distance “quantum fax machine” over optical fibers. 

Frequency correlated photon pairs can be generated in 
third order nonlinear waveguide by the spontaneous four 
wave mixing process (SFWM) [15,16]. Their state can be 
expressed as 

= j dnf{n)\uJp + n)s\uJp-^)^, (i) 

where ujp is the frequency of the pump light. The indices 
s and i indicate the signal photons (with frequency ujs) 
and idler photons (with frequency tOi), respectively. Q = 
uJs — ujp = ojp — LOi is the frequency detuning of signal 
photons or idler photons. /(Q) is the spectral amplitude 
of the biphoton state. 

The signal and idler photons are distributed to two 
parties, named Alice and Bob, over optical fibers. At 
Alice side, there is an object with a specific reflectivity 
pattern. The signal photons are firstly dispersed to dif¬ 
ferent directions along a line according to their frequen¬ 
cies by a spatial dispersion component, such as a grat¬ 
ing. Then they illuminate the object along a line, which 
is named as the illuminating line. Signal photons with 
different frequenciy will arrive at different positions on 
the object. The signal photons with a specific frequency 
detuning O would arrive at a specific position xyi in the 
line. The reflectivity pattern along the illuminating line 
would modulate the spectrum of the reflected signal pho¬ 
tons. Hence, the positive-frequency field operator of the 
signal photons at Alice side can be expressed as 

i?+(G,Ls)~ J dUlr{xn)as{ujp + ^) (2) 

where as{ujp+^) is the annihilation operator of the signal 
photons at the frequency ujp + Q. ts is the detection 



2 


time of the signal photons, Lg is the length of the optical 
fiber between the photon-pair source and the SPD at 
Alice side. Assuming that Lg is small, the group velocity 
dispersion (GVD) in the hber has been neglected. The 
phase coefficient of the signal photons in the fiber has 
been expanded in the vicinity of ujp + Vlo , which is the 
central frequency of signal photons, as j3g{ujp+fl) = /3s0 + 
— f^o)- 

At Bob side, the idler photons are temporally dispersed 
before detected by a SPD. The most convenient way to 
realize the dispersion is utilizing the GVD of the trans¬ 
mission fiber, by which the idler photons are sent to Bob 
side from the source. If the length of the optical fiber 
between the photon-pair source and the SPD at Bob side 
is Li, then the positive-field operator of the idler photons 
at time ti can be expressed as 

E+iU,L,) ~ J dna,{ojp - (3) 


where ai{uip — D) is the annihilation operator for idler 
photons at the frequency ujp — LI. The phase coefficient 
of the idler photons in the fiber can be expanded as 

Piiujp — LI) = /3io + Pii{Lio — D) -I- -l3i2{Llo — Ll)^, (4) 

where ujp — Oq is the central frequency of idler photons. 
The term with /3i2 is corresponding to GVD. Higher order 
dispersions in the fiber have been neglected. 

The signal and idler photons are detected by SPDs at 
both sides with arrival time recorded. Then Alice sends 
the single photon events to Bob. Time discriminated co¬ 
incidence measurement can be carried out at Bob side. It 
can be analyzed by the second-order Glauber correlation 
function Under the assumption that the temporal 

dispersion introduced at Bob side is large enough [see 
Section 1 of Supplement 1], the function can be ex¬ 
pressed as 


G^^\tg,Lg-U,L,) = |(0|A+(U,L,)4+(G,L,)|4')| (5) 


where r = (G - PgiLg) - (U - jduLi - fiaLloLi). 

According to equation (5), the coincidence measure¬ 
ment result has the profile of the spectrum of the bipho¬ 
ton state, which is modulated by the reflectivity pattern 
along the illuminating line in a nonlocal way [17,18]. 
Hence, the reflectivity pattern r{xQ) can be extracted 
from the coincidence measurement result. In this scheme, 
the signal photons are detected without discriminating 
their frequency and arrival time, hence, the reflectivity 
pattern along the illuminating line can not be recovered 
by Alice side only. On the other hand, the idler photons 
are discriminated according to their frequencies by the 
arrival time, supporting the nonlocal imaging process in 
the way of QGI. 

The principle can be explained by the transforma¬ 
tion of quantum correlation in this scheme. At Al¬ 
ice side, signal photons with different frequencies are 


mapped to different illuminating positions on the object 
due to spatial dispersion component. While, at Bob side, 
the idler photons with different frequencies have differ¬ 
ent arrival time due to the GVD of the hbers. Hence, 
the frequency correlation in the photon pairs is trans¬ 
formed to the correlation between the illuminating posi¬ 
tions of signal photons on the object and the arrival time 
of idler photons. Compared with QGI schemes based 
on momentum-momentum or position-position correla- 
tion[l] , our scheme realizes QGI in time domain thanks 
to the specihc correlation. Since the telecom-band fre¬ 
quency correlated photon pairs can be easily distributed 
over optical fiber, this scheme provides a practical way 
to realize long distance QGI. 



FIG. 1. (a). The frequency correlated photon-pair sonrce 

based on the SFWM process in a silicon nanowire waveguide. 
FI, F2 and F3 are optical filters for pump and signal/idler 
photons. EDFA: Erbium doped fiber amplifier, (b), The joint 
spectral density of the biphoton state. 


In the experiment, the frequency correlated photon 
pairs are generated by the SFWM process in a silicon 
nanowire waveguide [15]. The sketch of the source is 
shown in Fig. Ha)- The pulsed pump light has a cen¬ 
ter wavelength of 1550.92 nm, a pulse width of about 
3.7 ps and a repetitive rate of 40 MHz. It is gener¬ 
ated by a mode locked laser, then a filter (FI) is used to 
confine its spectrum and extend its pulse width. An er¬ 
bium doped fiber amplifier (EDFA) is used to increase the 
pump power and the filter (F2) after the EDFA is used 
to suppress the amplified spontaneous emission at the 
signal and idler wavelengths. The filters FI and F2 are 
composed of dense wavelength division multiplexing de¬ 
vices (DWDMs). The pump light is injected into a silicon 
nanowire waveguide (11 mm in length, with a cross sec¬ 
tion of 500 nm x 220 nm) through a lensed fiber. Thanks 
to the high nonlinearity and the designed low dispersion 
property of the silicon nanowire waveguide, broad-band 
frequency correlated photon-pairs can be generated. The 
signal and idler photons are separated by a filter system 
(F3) composed of cascaded coarse wavelength division 
multiplexing devices (CWDMs). The center wavelengths 
of signal and idler photons are 1530 nm and 1570 nm, 
respectively. 
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The joint spectrum density (JSD) of the biphoton state 
is measured by the method shown in ref. [19]. The signal 
and idler photons are dispersed temporally by a piece of 
50 km-long SMF. The measured JSD is shown in Fig. 
1(b). It is shown that the wavelengths of signal and idler 
photons are anti-correlated. The wavelength spans of the 
signal and idler photons are 16 nm, which are determined 
by the bandwidth of the filters F3. 


FIG. 2. Experiment setup. The frequency correlated photon- 
pair source (S) has been shown in Fig. [^The signal photons 
are distributed to Alice through a short piece of SMF and to 
Bob through 50 km-long SMF. At Alice side, the signal pho¬ 
tons are collimated to a spatial optical beam by a collimator 
(C) after an optical fiber circulator (CIR), then spatially dis¬ 
persed by a grating (G), and illuminate the object (O) after 
a focal lens (L). At Bob side, the GVD of the transmission 
fiber is utilized to realize temporal dispersion. At both sides, 
photons are detected by SPDs (SPDl & SPD2), and the time 
of single photon events are recorded(time recorder, TRl & 
TR2), and coincidence measurement (Goin.Meas.) is carried 
out. M: mirror. 

The sketch of the setup for long distance temporal QGI 
is shown in Fig. The source shown in Fig. [l] provides 
frequency correlated photon pairs. The signal photons 
are sent to Alice by a short piece of SMF. At Alice side, 
the signal photons are delivered to an spectrally encoded 
confocal microscopy (SECM) [20,21] system through a cir¬ 
culator (CIR). In the SECM configuration, the photons 
are collimated to a spatial optical beam with a 1/e^ di¬ 
ameter of 2.1 mm. A reflective diffraction grating (600 
Line/mm, blaze wavelength equal to 1600nm) is used to 
disperse the signal photons spatially. After a lens with a 
focal length I = 25.4 mm, the signal photons are focused 
on the surface of the object along a line, i.e., the illumi¬ 
nating line, with a length of about 250 pm . The signal 
photons with different frequencies would be focused on 
different positions along the illuminating line. 

The object is a standard photolithographic mask, 
which is a silica plate with a patterned chrome layer. 
It provides a reflectivity pattern. The signal photons are 
partially reflected according the reflectivity pattern along 
the illuminating line, resulting in modulated spectrum. 
The reflected signal photons are collected by the collima¬ 
tor and sent to a SPD (ID220, ID Quantique) through 
the CIR. The single-photon events are recorded with high 
time-resolution using a time recoder (TRl) in a photon 
correlator (DPC 230, Becker & Hickl CmbH). 

While, the idler photons are sent to Bob through 50 
km-long SMF. Then the idler photons are detected by the 


other SPD and the single photon events are also recorded 
with high time-resolution (TR2). On one hand, the fre¬ 
quencies of the idler photons are preserved after the prop¬ 
agation through such a long fiber, hence, the long dis¬ 
tance distribution of the frequency correlated biphoton 
states is realized. On the other hand, the transmission 
fiber provides a GVD of d = 900 ps/nm (estimated by 
the group velocity dispersion parameter of SMF). Since 
the bandwidth of idler photon is 16 nm, the idler photons 
are temporally broadened to 15.4 ns. 

When sufficient single photon events are collected at 
Alice and Bob sides, Alice should send the records of 
the detected signal photons to Bob by a classic channel, 
and Bob could realize coincidence measurement utilizing 
recorded signal and idler photon events. In the experi¬ 
ment, the single photon events at both sides are recorded 
by two channels of the photon correlator with a resolu¬ 
tion of 164.61 ps, and the coincidence is calculated by a 
computer according to the records. A typical measure¬ 
ment result is shown in Fig. [^a), while the inset shows 
the reflectivity pattern of the object. The dimension of 
the pattern is 150 pm x 140 pm. The red areas are the 
regions with high reflectivity. The dashed line is the il¬ 
luminating line corresponding to the histogram. It can 
be seen that the histogram has the shape of the pattern 
along the illuminating line clearly. 



FIG. 3. QGI of the object, (a), A typical coincidence his¬ 
togram. The inset is the reflectivity pattern of the object 
and the dashed line is the illuminating line corresponding to 
the coincidence histogram, (b). The image obtained by step- 
moving the object with a step of 10 pm. 

By step-moving the object along the direction orthog¬ 
onal to the illuminating line, two-dimensional imaging of 
the object can be obtained, which is shown in Figj^b). 
The moving step is 10 pm. The measurement time for 
each illuminating line is 10 minutes. The coincidence 
counts in every time bin are indicated by different colors. 
It can be seen that a clear image of the object has been 
obtained. Considering the long distance between Alice 
and Bob, this temporal QGI scheme realizes a function, 
of long distance image transmission like a “quantum fax¬ 
ing machine”. It is worth noting that in the experiment 
the image generated by the temporal QGI has a little 
distortion. It is due to the fiber length variation when 
the temperature changes during the measurement. Since 
the fiber is as long as 50 km, a temperature variation 
of 1°C would lead to an arrival time variation of several 
nanoseconds for the idler photons traveling through the 
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fiber. 

To analyze the resolution of our QGI scheme, another 
object with narrower line width is used, as shown in 
Figj^a). The image generated by the temporal QGI ex¬ 
periment is shown in FigQb). The line width of the 
object is 13 pm at the dashed line, while the width of 
the image at the corresponding position is about 1.65 ns. 
According to the experiment parameters [see Section 2 
of Supplement 1], it is related to a line width of 25.6 pm, 
much wider than the actual width. The broadening of 
the line is due to the limited resolution of the experiment, 
which is mainly determined by the wavelength resolving 
power of the grating at Alice side and the timing jitters of 
the single photon detecting system (including the SPDs 
and the photon correlator). 
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FIG. 4. QGI using a object with narrower line width, (a), 
The reflectivity pattern of the target object. The line width is 
13 um at the dashed line, while the width of the image at the 
corresponding postion is 25.6 pm. (b). The image generated 
by the temporal QGI experiment. 


Along the illuminating line, a reflectivity point at the 
object has a point-spread spectral function, which has a 
full width at half maximum (FWHM) determined by the 
resolving power of the grating [22], 

<5A = (6) 

where Aso = is the centeral wavelength of signal 

photons in vaccum, N is the total number of grooves in 
the beam area on the grating. The point spread spectral 
function will result in a point-spread coincidence peak, 
with a FWHM of 


6 ti = 6 X X d, (7) 

where d is the temporal dispersion at Bob side. On the 
other hand, the impact of the single photon detecting 
system is denoted by St 2 , including the effect of timing 


jitters of two SPDs and the time resolution of photon 
correlator. Hence, the resolution of the temporal QGI 
can be estimated by 


a/ + 6 t^ d 0 

= - i -'dX 




p COS 9( 




where, 0 o is the diffraction angle of the signal photons 
with wavelength Aso, P is the period of the grating. The 
relation p(sin0 —sin = A between the diffraction angle 
0 of the signal photons, the wavelength A and the incident 
angle 0^, has been used in the equation. 

According to the parameters of the experiment setup 
{I = 0.0254 m, p = 1.67 pm, 0o = 11.9 , A^o = 1530 nm, 
N = 1044 grooves, 6 t 2 = 389 ps, d = 900 ps/nm), the 
calculated spatial resolution is about 23.8 pm, agrees well 
with the experiment results considering the fact that the 
reflectivity of the object has a rectangluar profile. 

Equation (7) shows that the spatial resolution could 
be improved through (1) increasing temporal dispersion 
at Bob side; (2) improving the time resolution of the 
single photon detection system; (3) increasing the beam 
diameter at Alice side; (4) utilizing an objective lens with 
smaller focal length at Alice side. Our recent work [23] 
shows that T 2 = 80 ps can be realized in the single photon 
detection system composed of super conducting nanowire 
single photon detectors (timing jitter ^ 60 ps) and high 
performance photon correlator (timing jitter ~ 12 ps). If 
the diameter of the signal photon beam is expended to 5 
mm and an objective lens with a focal length of 5 mm is 
used, a resolution as high as 1.1 pm can be expected. 

As a conclusion, we proposed and experimentally 
demonstrated a scheme of temporal QGI. By step-move 
the object, the object can be imaged through 50 km SMF 
by this scheme. This scheme works like a quantum fax 
machine over optical fibers. The spatial resolution of this 
scheme is analyzed, showing that it could be less than the 
wavelength of the illuminating photons by proper system 
design. It will extend the application of QGI at large ge¬ 
ographical scale. 
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^Appendix 


SI. function 


It is well known that for biphoton state Idt), the second- 
order Glauber correlation function can be calculated 
using [1] 




{0\Et{ts,Ls)E+{U,L,m 


,( 9 ) 


where tg, ti is the detection time of signal and ider pho¬ 
tons, La^ Li is the optical path of signal and idler pho¬ 
tons, respectively. Including the dispersion of the trans¬ 
mission fiber, phases will be introduced to the positive- 
frequency operators, as shown in the Eq.(2) and Eq.(3), 
following which {Q\E'^{tg, Ls)Ef {U, Li)\'^) can be calcu¬ 


lated with the state in Eq.(l), 

{Q\Et{ts.Lg)Ei{ti,L,m 

X (0| 0-5 \ ^p H“ (0 (^i) \ ^p i 

xS{u!p + fl — UJg)d{ujp — rt — LOi) (II) 

= J ( 12 ) 

= 27re^‘^EAfmrixn)) * (13) 


where we used the phase coefficient expandation in 
the Eq.(2) and Eq.(4) and the commutation relation 
[UmiuJl), ai{uj 2 )] = 6m,nS{uJl - UJ 2 ), rn,n & {s, i}, Sm,n is 
the Kronecker delta function, and S(-) is the Dirac delta 
function, includes all the D-independent phases, and 
T = {tg- PgiLg) - {ti - PiiLi - f3i2LloL^). Fr{-) denotes 
the inverse Fourier transformation, as 


Er{.f{^)r{xn)) = ^ y dD/(U)r(a;n)e^^'^, 

jr^(g-ift2nV2) = — f 

27r J 


(14) 


(15) 


Combining equation (13), (14) and (15), we obtain 


{Q\Et{tg,Lg)Et{ti,Li)\^) 

- y drily dD/(D)r(a;n)e^'^^y 
~ y dTi|y dD/(D)r(a;o)e'^^'^y (17) 

~/(0)r(a;n)|n=^- (18) 

From Eq.(17) to Eq.(S18), we have assumed that a very 
large dispersion has been introduced at Bob side, so that 


PaLi is much larger than the temporal width of idler- 

^2 

photon wavepackets, and neglected the term , be¬ 

cause 2 ;^ < 1 [2] • 

Lingking Eq.(9) and Eq.(18) for clearity, we can get 
the expression of G^^i function when large temporal dis¬ 
persion is introduced at Bob side, 

G^ \tg,Lg]ti,Li) ~ l/(ll)^(2;o)|f2=T//3i2ii > (^9) 

where t = {tg - PgiLg) - {U - PnL^ - PaLl^Li). 

According to Eq.(19), on one hand, it is indicated that 
the coincidence measurement results have the shape of 
the spectrum of the biphoton state, which is nonlocally 
modulated by the reflectivity pattern of the object, so by 
extracting the reflectivity spectra r{xn), the image of the 
object can be reconstructed. On the other hand, it can be 
seen that, after temporal dispersion, the coincidence peak 
is spread to have a width of r^, = PaLiLlui-, where 0.^ is 
the spectrum width of the biphoton state. The spreading 
of the coincidence peak is due to the GVD introduced at 
Bob side. For idler photons, using the expandation of 
the phase coefficient in Eq.(4), the group velocity of idler 
photons can be expressed as [2] 

— = = (20) 

Vig OUJi OS I 

So the time needed for idler photons from the source to 
the SPD at Bob side will be 

tn = — = PiiLi + Pi2Li{VlQ — n). ( 21 ) 

Vig 

For idler photons with frequency difference AD, the time 
delay difference will be 

At = -Pi2L,An, (22) 

which is coincident with the spread of the coincidence 
peak. In fact, after the spatial dispersion at Alice side, 
and the temporal dispersion at Bob side, the frequency 
correlation in the photon pairs is transformed to the cor¬ 
relation between the illuminating position X(i of signal 
photons on the object and the travel time tQ of idler 
photons before detected by the SPD. Compared with pre¬ 
vious QGI schemes based on momentum-momentum or 
position-position correlation, our time domain QGI is re¬ 
alized utilizing this kind of correlation. 


S2. SCALE TRANSFORMATION RELATION 
BETWEEN THE OBJECT AND THE 
TIME-DOMAIN IMAGE 

This section discusses the relation between the scale of 
the object and the image. It is clear that in the direction 
orthogonal to the illumiating line, the scale of the image 
should be equal to the scale of the object. In the follow¬ 
ing, we focus on the scale transformation relation in the 
direction of the illuminating line. 
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At Alice side, the relation between the illuminating 
position XQ on the object and the signal photon frequency 
Wp + n can be expressed as 


_ - Oo d6» 

(c. + nofdA 


(23) 


\—2ttcI (cjp+Qo) 


where xho is the illuminating position of signal photons 
with frequency a;p+Oo(the central frequency of the signal 
photons), I is the focal length of the lens, c is the velocity 
of light in vacuum, 0 is the diffraction angle of light beams 
after grating, and A is the wavelength of signal photons, 
which is corresponding to the frequency by the relation 
A = 27rc/(wp+n). As we use the first order of interference 
of the grating, the angle of the incident light 9i, the angle 
6 of the diffration light with wavelength A, must satisfy 
the equation [3], p (sin0 — sin6*i) = A, where p is the 
period of the grating. Hence, 


dA 


1 


A=- 




p COS 00 ’ 


(24) 


where, 0o is the diffraction angle of signal photons with 
frequency Wp + Hq. Two points with high reflectivity on 
the object with distance Aa; will reflect signal photons 


with frequency differency equal to 

AH = — ^o) ^ 0qAx. (25) 

where, 0q is the diffraction angle of signal photons with 
frequency Wp + Hq . 

According Eq.(22), the time differency between the two 
coincidence peaks resulted from the two high-reflectivity 
points will be 


At = -PaLi/AQ, = ) ygcos 6 * 0 Ax. (26) 

Zircl 

Considering the relation between the dispersion pa¬ 
rameter of d and Pi 2 , i.e. d = —A^o = 

''iO 

2ncl{ujp — Hq) is the central wavelength of idler photons, 
equation (26) can be expressed using the dispersion pa¬ 
rameter d as 


A2 

At = —^^pd cos 0o^x, (27) 

where A^o is the central wavelength of signal photons, 
and Aso = Zucj{ujp + Hq). 
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